Interaction of a Contact Resonance of Microspheres with Surface Acoustic Waves 
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We study the interaction of surface acoustic waves (SAWs) with the contact-based, axial vibra- 
tional resonance of 1 /xm silica microspheres forming a two-dimensional granular crystal adhered 
to a substrate. The laser-induced transient grating technique is used to excite SAWs and measure 
their dispersion. The measured dispersion curves exhibit "avoided crossing" behavior due to the 
hybridization of the SAWs with the microsphere resonance. We compare the measured dispersion 
curves with those predicted by our analytical model, and find excellent agreement. The approach 
presented can be used to study the contact mechanics and adhesion of micro- and nanoparticles as 
well as the dynamics of microscale granular crystals. 

PACS numbers: 68.35.Iv, 78.47.jj, 46.55.-fd, 45.70.-n 
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Wave phenomena in granular media is a rich and 
rapidly developing field of research . At the heart of 
this field is the Hertzian model of elastic contact between 
spherical particles, in which the stiffness of the contact 
depends on the applied force One type of granular 
media, often referred to as "granular crystals", consists 
of close-packed, ordered arrays of elastic particles that 
interact via Hertzian contact 0, IH. Granular crystals 
have been shown to support a wide range of linear and 
nonlinear dynamical phenomena not encountered in con- 
ventional materials, and have been suggested for various 
engineering applications 0, [3| • 

Acoustic studies of granular media typically in- 
volve macroscopic particles with dimensions of ^ 0.1- 
10 mm [iHl]; whereas contact-based vibrations of mi- 
croparticles with dimensions of (or under) ~ 1 /xm re- 
main largely unexplored. The scale factor is significant 
as a microparticle system cannot be thought of simply 
as a scaled down version of a macroscale system which is 
governed by the same physics. Rather, microparticles are 
expected to yield qualitatively different dynamics. One 
crucial factor is the role of adhesion 0, Q , which is almost 
negligible on millimeter scales but significant on micron 
scales. Because of adhesion, a microsphere in contact 
with a substrate is pulled toward the latter. This results 
in an equilibrium contact stiffness and an axial "contact 
resonance" vibrational mode with frequency determined 
by the particle mass, the adhesion, and the elastic prop- 
erties of the particle and substrate [1] . This phenomenon 
has not hitherto been observed experimentally, although 
rocking-mode vibrations at much lower frequencies have 
been studied in ~ 5-50 ^m spheres [3-3] , and vibrations 
corresponding to free particle eigenmodes at much hi ghe r 
frequencies have been studied in 5 yitm silica spheres |10| . 

In this letter, we study the contact-based, axial vibra- 
tional resonance of 1 /im diameter microspheres forming 
a two-dimensional granular crystal adhered to a substrate 
using another hitherto unexplored phenomenon, i.e., the 
interaction of axial contact resonances of microparticles 
with surface acoustic waves (SAWs) in the substrate. 



We use the laser-induced transient grating (TG) tech- 
nique (111, II4I to excite the long-wavelength (relative to 
the particle size) SAWs and measure their dispersion. 
The measured dispersion curves exhibit classic "avoided 
crossing" behavior due to the hybridization of the SAWs 
with the contact resonance of microspheres. Such cou- 
pling between SAWs and mechanical surface oscillators 
was studied in theoretical works [13l - ll5{ , and the experi- 
mental work [l6'|. We analyze our measurements using a 
simple model that yields an analytical expression for the 
dispersion relation. We find excellent agreement between 
our model and the measurements using the frequency 
of the contact resonance as a single fitting parameter, 
and we compare our results with estimates based on the 
Derjaguin-MuUer-Toporov (DMT) contact model [13] • 
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FIG. 1: [Color online] (a) Microspheres interacting with a 
SAW via contact "springs", with notations used in the theo- 
retical model, (b) Photograph of the sample, (c) Representa- 
tive image of the silica microsphere monolayer, (d) Schematic 
illustration of the TG setup. 



Our sample is a hexagonally packed monolayer of 
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D = 1.08 /im diameter silica microsplieres deposited on 
an aluminum-coated fused silica substrate, as shown in 
Fig. [TJb,c). The fused silica slab is 1.5 mm thick, and 
the aluminum layer, which serves as a medium to absorb 
pump laser light, is 0.20 /im thick. To assemble the mono- 
layer on the substrate, we used the "wedge-shaped cell" 
convective self-assembly technique l^- 2^ ■ The resulting 
monolayer has an overall area of ~ 5 mm by 10 mm, as 
shown in Fig.[ljb). A representative optical image of the 
monolayer packing is shown in Fig. [TJc). Although the 
monolayer has a relatively uniform distribution, defects 
and grains of uniform packing are clearly present. 

We use a laser-induced TG technique ll|, 12 1 to mea- 
sure the phase velocity dispersion of SAWs in our sam- 
ple. The TG setup with heterodyne detection used for 
these experiments has been described previously 12]. In 
summary, two excitation beams derived from the same 
laser source (515 nm wavelength, 60 ps pulse duration, 
2.44 /iJ total pulse energy at the sample) enter the sam- 
ple through the transparent silica substrate, as shown 
in Fig. [lid), and are overlapped at the aluminum layer 
forming a spatially periodic interference pattern. The 
pump spot has 500 /zm diameter at 1/e^ intensity level. 
Absorption of the laser light by the aluminum film in- 
duces rapid thermal expansion, which leads to the gener- 
ation of counter-propagating SAWs with wavelength As 
defined by the period of the interference pattern. The pe- 
riod is controlled by a phase mask pattern used to create 
the two excitation beams, by splitting the incident beam 
into +/ — 1 diffraction orders [uj. Switching phase mask 
patterns allows measurements at multiple acoustic wave- 
lengths. 

SAW detection is accomplished using a quasi-cw probe 
beam (532 nm wavelength, 10.7 mW average power at the 
sample) focused at the center of the excitation pattern to 
a spot of 150 fiva. diameter. The probe beam also enters 
the sample through the silica substrate and is diffracted 
by surface ripples and refractive index variations in the 
substrate induced by SAWs [2l|. The diffracted beam is 
overlapped with the reflected reference beam (local oscil- 
lator) [ill , [l^ and directed onto a fast avalanche photo- 
diode with a 1 GHz bandwidth. The signal is recorded 
using an oscilloscope and averaged over 10"* repetitions. 

Figure [2Ia,b) shows typical signal waveforms acquired 
at an acoustic wave vector magnitude k = 2ti/\s = 
0.46 /im~^. Figure [21 a) shows the signal from a sam- 
ple location without microspheres, and Fig. [^b) corre- 
sponds to a location with spheres. In both cases there is 
a sharp initial increase, which corresponds to the excita- 
tion pulse arriving at the sample. The slowly decaying 
component is due to the "thermal gratin g" a ssociated 
with the temperature profile in the sample [111 IT^. The 
high frequency oscillations are due to acoustic waves. 

Figure [2lc,d) shows the Fourier spectra of acoustic os- 
cillations corresponding to the signals in Fig. [2Ia,b) (2^ . 
In the off-spheres case there are two clear peaks, which 
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FIG. 2: [Color online] Normalized signal (a) off-spheres, and 
(b) on-spheres, for k — 0.46 fim~^ . The acquired signal S 
is normalized by the maximum signal amplitude Sq. Fourier 
transform (FT) magnitudes corresponding to the signals in 
(a) and (b), plotted in log scale (c) and linear scale (d). The 
black curve corresponds to the signal in (a), and the red curve 
to the signal in (b). The markers denote the identified peaks, 
which are plotted in Fig. [3] using the same markers. 



correspond to a Rayleigh SAW (the low frequency peak) 
and a longitudinal wave in the substrate 22] . Figure [31 
shows the acoustic dispersion curves. For the off-spheres 
case, we see linear dispersion curves that agree very well 
with lines corresponding to the longitudinal and Rayleigh 
wave velocities in fused silica. We used typical wave 
speeds for fused silica of cl — 5968 m/s (longitudinal) 
and ct = 3764 m/s (transverse) [23|, and we calculated 
the Rayleigh wave speed cr = 3409 m/s by numerically 
solving the Rayleigh equation [2J] . More accurate calcu- 
lations accounting for the aluminum layer pol [23 | showed 
that the reduction in the Rayleigh wave speed due to the 
aluminum layer does not exceed 1.4%. 

The on-spheres case yields starkly different behavior 
from the off-spheres case, as can be seen by comparing 
the signal waveforms shown in Fig. 2(a,b). The compar- 
ison of the spectra in Fig. 2(c) shows that the longitu- 
dinal peak is unaffected by the presence of the spheres 
whereas the SAW peak is split in two. The on-spheres 
dispersion curves in Fig.[3]reveal a classic "avoided cross- 
ing" between the Rayleigh wave and a local microsphere 
resonance. The lower branch starts as a Rayleigh wave at 
low wave vector magnitudes and approaches a horizon- 
tal asymptote corresponding to the resonance frequency. 
The upper branch is close to the Rayleigh wave at high 
wave vector magnitudes; in the "avoided crossing" region 
it deviates from the Rayleigh line and stops at the thresh- 
old corresponding to the transverse acoustic velocity of 
the substrate. 

To describe SAW propagation in our system of mi- 
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FIG. 3: [Color online] Dispersion relations. Red and black 
markers are the measured frequency peaks for the on- and off- 
spheres cases, respectivly. The solid red line is the dispersion 
calculated using our model. Also shown are lines correspond- 
ing to longitudinal, transverse, and Rayleigh waves in fused 
silica, and a horizontal line corresponding to the microsphere 
contact resonance frequency. 



crospheres coupled to an elastic substrate, we adapt an 
approach developed in the theoretical works l3-15|. As 
shown schematically in Fig.[IJa), we model the substrate 
as an elastic half-space (z < 0), where the surface of the 
halfspace is connected to an array of linear surface os- 
cillators, with mass m and stiffness K^-, which represent 
the microspheres connected to the substrate via "con- 



tact springs" . We calculate the mass of a microsphere as 
TO = |7r(-j)^ps, where ps = 2.0 g/cm^ is the density of 
the silica microspheres as provided by the manufacturer 
(Corpuscular, Inc.). 

We approximate the microspheres as point-masses, as 
the lowest spheroidal resonance of the microspheres, /2 = 
2.9 GHz [25|, is much greater than the frequencies ob- 
served in the experiment. The equation of motion for the 
surface oscillator can be written as mZ + K2{Z ~ Uz,o) = 
0, where Mz.q is the displacement of the substrate sur- 
face, and Z is the displacement of the oscillator relative 
to the surface. The particles exert a vertical force on the 
substrate, leading to the following boundary conditions 
at the surface z = 0: 



K2{Z-u 



z.O) 



A 



Oxz = 0, 



(1) 



where cr^^ and a^z are components of the elastic stress 
tensor |23] , and A — is the area of a primitive unit 

cell in our hexagonally packed monolayer. 

Since the acoustic wavelength is much larger than the 
sphere size, in Eq.[T]we use an effective medium approach, 
and approximate the average normal stress at the surface 
as the force exerted by the microsphere "spring" divided 
by the area of a unit cell (assuming perfect packing) even 
though the stress is localized around the contact area. 
We follow the standard procedure for the derivation of 
the Rayleigh wave equation Q , but use Eq. [1] instead 
of stress-free boundary conditions to obtain the following 
dispersion relation for the SAWs in our coupled oscillator 
system [i^: 
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where wo = 27r/o = ^jK^jTn, is the angular frequency of 
the contact resonance, and p2 = 2.2 g/cm'^ f26| is the 
density of the silica substrate. On the left-hand side of 
Eq. [2 the term in square brackets is familiar from the 
Rayleigh equation, and the term in the parentheses de- 
scribes the resonance of the oscillators. The right-hand 
side of Eq. [2] is responsible for the coupling between the 
Rayleigh waves and the oscillators; if it is made to vanish 
(for instance by assuming a vanishing areal density mj A)^ 
then the oscillators and SAWs in the substrate are effec- 
tively decoupled. We also see that A and p2 relate to 
the coupling strength, K2 relates to the frequency of the 
avoided crossing through wq, and to relates to both. 

By taking the frequency of the contact resonance as a 
fitting parameter, and using least squares minimization 



between the numerical solution of Eq.[2]and the measured 
dispersion, we find /o = 215 MHz. The fitted resonant 
frequency is plotted as the blue dashed line in Fig. [H 
and gives a contact stiffness of K2 = 2.7 kN/m. We plot 
the numerical dispersion curve from our fitting as the red 
solid line in Fig. [3] 

Any real solution of Eq. [5] must yield the phase ve- 
locity w/fc smaller than ct, otherwise at least one of 
the square root terms becomes imaginary. Therefore, 
the calculated upper dispersion branch terminates at the 
threshold w = ct^, in agreement with the experiment. 
In some cases leaky wave solutions with complex w can 
be found above the threshold , however we did not in- 
vestigate complex solutions since in our experiment the 
upper branch peak disappeared past the threshold. 
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Using the oscihator equation of motion, we estimated 
the relative displacements (and phase) of the micro- 
spheres and the surface for various limiting cases [20| . 
as is shown in Fig. [31 In the flat dispersion region of the 
lower branch, there is predominantly sphere oscillation 
with very small surface displacements. Indeed, our exper- 
imental data show that the amplitude of the correspond- 
ing peak becomes progressively smaller compared to the 
Rayleigh-like upper branch as the wave vector magnitude 
is increased [20|. 

We estimate the frequency of the contact resonance 
(and the contact stiffness) of the microspheres using the 
DMT contact model , and compare with the res- 

onant frequency obtained from the measured data. The 
DMT model describes the contact between an elastic 
sphere and a flat substrate under the presence of adhesive 
forces. The model assumes that the deformation profile 
is Hertzian and the adhesive forces act outside the con- 
tact area. For the case of small displacements, the full 
DMT model can be approximated as 17l |: 



F = KR^/^a^/^ - 2TrwR, 



(3) 



where _F is a force applied to the sphere, a is the dis- 
placement of the center of the sphere towards the sub- 
strate, w = 0.094 J/m^ is the work of adhesion between 
silica and alumina (as the aluminum surface is normally 
oxidized) 20], R is the effective microsphere radius of 
contact (which reduces to i? = Z3/2 for a sphere on a flat 

substrate [4]), and K = (|(^^ + ^^))"^ is the the 
effective modulus, where the aluminum has elastic mod- 
ulus El = 62 GPa and Poisson's ratio i^i = 0.24 [13], 
and the microspheres have elastic modulus Eg = 73 GPa 
and Poisson's ratio — 0.17 [i^. Using Eq. [3l we calcu- 
late the equilibrium displacement uq — ( 27ru^_^-j2/3 _ 
0.44 nm, the linearized stiffness around the equilibrium 
point K2,DMT = f (27rwi?2if2)i/3 ^ 1,1 kN/m, and 

^^/K2,DMT/m 



140 MHz. Below the ax- 



fo,DMT - 2^ 

ial contact resonance, a rocking mode is also predicted 
frock = WT^^I^ = 10 MHz @, however this is sig- 
nificantly below our measured acoustic frequency range. 
The discrepancy between the estimated and the mea- 
sured values of /o may be caused by uncertainties in the 
contact and adhesion models. Challenges in application 
of DMT to real nanoscale contacts have been underscored 
by studies in atomic force microscopy 0, [2^ . Typical ad- 
hesion studies relying on measuring a pull-off force pro- 
vide limited information for verification of adhesion the- 
ories 0, Q . Our experiment provides a direct pathway to 
the contact stiffness, and thus offers a promising tool for 
studying adhesion and contact mechanics. 

Potentially useful information can also be obtained by 
studying the attenuation of the SAWs. As can be seen in 
Fig. [2][d), the peaks in the on-spheres case are broader. 
Faster attenuation is also evident from a comparison of 
the signal waveforms in Fig. [Ha,b). In the off-spheres 



case, the acoustic signal decays as the SAW wavepacket 
leaves the probe spot, with material losses bein g ii egligi- 
ble within the time window of the measurement [30| . The 
apparent decay time is thus determined by the width of 
the wavepacket (i.e. the width of the excitation spot) 
and the SAW group velocity. In the on-spheres case the 
group velocity is lower, yet the decay time is shorter, 
which indicates that the presence of the spheres leads 
to additional attenuation. One possible mechanism is 
scattering due to the spheres. A single oscillator on the 
surface of a halfspace will radiate acoustic energy into 
the substrate; however, a collective mode of a periodic 
array with w < c^/e will not radiate 15, 3^. In our case, 
the sphere packing is not perfectly periodic, which may 
lead to scattering and radiation into the bulk. In the flat 
region of the lower dispersion branch, the acoustic mode 
is close to the contact resonance of the spheres. In this 
case, the peak width may also be determined by inhomo- 
geneous broadening caused by sphere-to-sphere variation 
of the contact stiffness. If this is the case, then the peak 
width can be used to study the statistics of contact stiff- 
nesses. Peak broadening may also be caused by anhar- 
monicity, but we estimate sphere displacements to be in 
the linear regime [l^ , and no anharmonic effects such as 
second harmonic generation are observed. In the future, 
it would be interesting to study the nonlinear regime, 
which should be achievable based on observations of a 
complete detachment of micron-sized particles by high- 
amplitude SAWs [31]- 

We have seen that an avoided crossing between the 
Rayleigh wave and contact resonance of the spheres oc- 
curs at wavelengths much larger than the sphere size and 
is well described by the effective medium approximation. 
Thus our structure belongs to a class of "locally resonant 
metamaterials" , for which interesting effects have been 
observed on the macroscale 32, 3^1 ■ We expect the effec- 
tive medium model to break down at shorter wavelengths 
where phononic crystal effects should be seen [3J]. Fur- 
thermore, our model treats the spheres as independent 
oscillators that interact through the elastic substrate but 
not directly. Although the spheres are closed-packed, the 
model describes the data very well. We believe that this 
is, again, due to the fact that the acoustic wavelength 
is large compared to the sphere size. At shorter wave- 
lengths we expect to see rich dynamics due to interaction 
between the spheres 35, 3^; the interaction with SAWs 
will make such dynamics even richer. 

In summary, we have studied the interaction of SAWs 
with the contact-based resonance of microspheres form- 
ing a two-dimensional granular crystal. The experimen- 
tal method can be used to study the adhesion and contact 
mechanics of microparticles. It also enables the study 
of granular crystals on the microscale. A rich array of 
dynamical phenomena observed in macroscale granular 
crystals, and their promise for practical applications 
suggest interesting possibilities for microscale granular 
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crystals. An analogy can also be made between SAWs in- 
teracting with a microsphere contact resonance and sur- 
face plasmon-polariton waves in a metal film interacting 
with a localized surface plasmon resonance of a metallic 
nanoparticle coupled to the metal film 37[. This may 



lead to acoustic analogies of some plasmonic phenomena 
and applications (38| . Finally, the nonlinearity of the 
Hertzian contact holds promise for an application of our 
approach to developing nonlinear SAW devices. 
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N.X.F. are also partially supported by NSF grant CMMI- 
1120724. 
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SPECTRA FOR MULTIPLE SAW WAVELENGTHS 

In Fig. [U we show the measured spectra, as in Fig. 2(c) of the main text, for multiple SAW wavelengths. We 
measure each signal at a 20 GHz sampling rate, and normalize by the maximum amplitude of the signal. Starting 
after the sharp initial increase, the slowly decaying component is identified using a smoothing function and subtracted 
from the normalized signal. Using this signal with the slowly decaying component subtracted, we select a window of 
0.24 /is duration. This window was then zero-padded to yield a window of 0.36 /iS total duration, which was used to 
calculate the FT magnitudes. In some cases, when a small peak is located on the "shoulder" of a larger peak, the 
appearance of the small peak depends of the exact choice of the FT time window: depending on the starting position 
of the window, the small peak may appear either as a positive peak or a "dip" . Accordingly, we adjusted the start 
time (by up to 15 ns after the sharp initial increase, for each wavelength) to ensure that small peaks appear as positive 
ones. 
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FIG. 1: Measured spectra for multiple SAW wavelengths. 
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SAMPLE FABRICATION DETAILS 



The fused silica wafer was cleaned using Piranha solution prior to aluminum deposition. Following the cleaning, 
a 0.20 /im thick aluminum coating was deposited via electron beam evaporation. The aluminum-coated substrate 
was then spin-coated with 52-14 photoresist, and diced into 2 cm square pieces. The photoresist was removed by 
submerging the substrate in acetone and then isopropanol baths, and dried under N2 flow. 



To assemble the monolayer on the surface, we used a "wedge-shaped cell" convective self assembly technique |lj, [2j. 
We used 1.08 /im diameter silica microspheres, purchased from Corpuscular Inc. as a suspension of 5.0 wt% in water. 
Prior to the deposition, the suspension was further diluted to 1.25 wt%. Following [J, we submerged the substrate 
and two glass microscope slides in a 30% hydrogen peroxide bath, at 80° C for 1 hour. The substrates and slides were 
then cleaned with acetone, isopropanol, deionized water, and dried under N2 flow. The substrate was then epoxied 
on top of one of the glass slides, at a distance from the edge such that, when the second glass slide was clamped on 
top of this assembly at one edge, the angle of the wedge formed between the substrate and the top glass slide is 3 
degrees. We then deposited 20 fiL of the diluted suspension on top of the aluminum-coated substrate, and clamped 
the second glass slide on top (by the edge) to form the wedge. The whole setup was then placed on a tilt, such that 
the meniscus drying front will recced upwards on an incline of 10 degrees. 



Following [3|, we calculated the wave speed for a generalized Rayleigh wave in an elastic halfspace (silica) coated 
with a thin finite layer (0.20 fiui thick aluminum) on the top surface. The solution for the layered halfspace allows for 
higher order surface waves, but for our range of SAW wavelengths we are below the cutoff point for the existence of 
these additional solutions. Additionally, in this case, the Rayleigh wave speed is dependent on the SAW wavelength. 
We plot this relationship in Fig. [21 for our substrate parameters. At the smallest wavelength (highest difference in 
wave speed from the pure fused silica halfspace), the wave speed is a maximum of 1.4% slower than in pure fused 
silica. 




RAYLEIGH WAVE SPEED IN LAYERED SUBSTRATE 



1 




0.98 



0.2 



0.4 0.6 
wavenumber |iim~^] 



0.8 



FIG. 2: Rayleigh wave speed cr^l in a fused silica substrate coated with 0.20 fim thick aluminum film, normalized by the 
Rayleigh wave speed in pure fused silica cr, for the experimentally measured wavenumbers. 
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DERIVATION OF SAW DISPERSION RELATION FOR AN ELASTIC HALFSPACE WITH COUPLED 

SURFACE OSCILLATORS 

Elastic waves in the substrate are described by wave equations V'^cj) = '^^W ^^"^ V^V = 73- ^^T: where t is time, 
and 4> and ip are dilational and transverse potentials related to the displacement components and Uz as follows: 



M.7 



— — — — and Uz = §7 + §7 [sl- Assuming a traveling wave solution of the form (p — /(z)e''^'^* '^^^ and 



dz ^ dz dx 



■0 = g{z)e^'^'^*' we obtain solutions to the wave equations of the form 



kz. 1 y j -^i(iut — kx) 

Bie V " "L 

(1) 



kz. 1 ^—r^-\'i(ujt — kx) 

The equation of motion for the surface oscillator can be written as: 

mZ + K2iZ ~ Uz,o) ^ 0. (2) 

We assume a traveling wave solution for the oscillator motion of the form Z = e'^"*^*^^), which we substitute Eq. [2] 
to obtain: 

Z^4^- (3) 
Substituting Eq. [3]into the boundary conditions at z = (shown in Eq. 1 of the main text), we obtain: 

Using Eq. [l] Eq. |4j and isotropic linear elastic stress-strain relations , we obtain the dispersion relation for the 
SAWs in our coupled oscillator system, as shown in Eq. 2 of the main text. Additionally, we note that Eq. [3] was 
also used to estimate the relative oscillator and surface displacements (and their relative phase), as shown in Fig. 3 
of the main text. 



WORK OF ADHESION AT A SILICA-ALUMINA INTERFACE 



We estimate the work of adhesion w between our silica microspheres and the aluminum-coated silica substrate. As 
the aluminum-coated substrate has a thin oxidized layer (alumina) on top of the aluminum, we calculate the work of 
adhesion at a silica-alumina interface. Following 4] , we find the work of adhesion in terms of the Hamaker constant 
A12 for a silica-alumina interface, where w = i2nl)i (Dq = 0.165 nm is a standard value used for the interfacial cutoff 
separation distance for a variety of media Q). This coefficient takes into account the van der Waals forces between 
the two surfaces in contact. The Hamaker constant can be defined in terms of the material properties of the two 
materials in contact, interacting across a third medium (in this case air), such that: 

^ 3fcr £1-63 £2-e3 ^ 3/it/e ("1 - nDjnl - n|) 

4 ^ei+e3^^e2 + e3^ 8^/2 (n2+n§)i/2(n2+n2)i/2[(n2 + „2)i/2 + („2 + „2)i/2]' 

where subscripts 1, 2, and 3 correspond to silica, alumina, and air, respectively, T = 293 K is room temperature, 
Ve = 3.2 • 10"'^^ Hz is the main electronic absorption frequency in the UV, k is Boltzmann's constant, h is Planck's 
constant, permittivities ei = 3.8 F/m, 62 — 11 F/m, and £3 = 1.0 F/m, and refractive indices rii = 1.45, n2 = 1.75, 
and n3 = 1.00. All of the preceding constants are tabulated in 4^. Using Eq. [Sj we calculate A12 = 9.6 • 10"^" J. This 
gives a work of adhesion oi w = 0.094 J/m^. 



CONTACT MODEL 



There are several models that are commonly used to describe the contact of elastic spheres under the presence of 
adhesive forces This includes the DMT [6], Jordan-Kendall-Roberts (JKR) [7], and Mau gis-Dugdale [8| models. 
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In the main text, we used the DMT model to estimate the frequency of the contact resonance, which assumes that 
the deformation profile is Hertzian and the adhesive forces act outside the contact area @. In the JKR model, the 
deformation profile is not constrained to be Hertzian, and the adhesive forces are assumed to act only within the 
contact area Q. The DMT model is typically applied to weakly adhesive systems with small, stiff particles, and 
the JKR model to strongly adhesive systems with large, soft particles Q. The moregeneral Maugis-Dudgale model 
continously transitions between these two regimes, and can be written as follows 0, 0|: 



1 = ^)2/3[^M2 - 1 + (A/2 - 2)(arctanv/M2 - 1)] 



4X^a, K 



(— -M+ \Jm'^ - l(arctan\/M2 - 1)1, 



F 



3 nR'^w 



Ck = — — 



_ Xa^ ( ) 1/3 [T^giTT + (arctan VM^)] , 



R ' R 



where a is the radius of contact, AI is the width of an annular region where adhesive forces are assumed to act, 
zq = 0.4 nm is the approximate interatomic distance at the interface and A = ^^(^^)^^'^ ^ non-dimensional 
parameter that describes the behavior of the contact. Equation [S] approaches the JKR model as A — oo, and the DMT 
model as A — 0. For our silica microspheres in contact with an aluminum-coated substrate, we calculate A = 0.5. We 
estimate the frequency of the contact resonance using Eq.[6l and find it to be less than 5% smaller than the frequency 
predicted using the DMT contact model. 



SURFACE DISPLACEMENT AMPLITUDE 

Following we estimate the displacement amplitude in the laser-generated SAW to be — g^Qo = 23 pm, 
where: 

_ 3aLil~ Rref) '^cl^ 

r = [(4 - - + (ci - 4)-i] - 2]-\ (8) 

Qo — QlHt^R^) is the laser energy density (with a pulse energy of = 2.44 ^J and spot radius of Rl — 150 /im), 
and = 23.6 /im/m-K is the linear thermal expansion coefficient, pi — 2.7 g/cvp? is the density, Cp = 0.90 kJ/kg-K 
is the specific heat, and Rref = 0.92 is the surface refiectivity of aluminum [10]. Although the microsphere vibration 
amplitude can be greater than the substrate surface displacement, this value is over an order of magnitude less than 
the equilibrium displacement of the microsphere due to adhesive forces, and thus serves as a basis for assuming the 
linearity of the response. 
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